Effects of competition on pattern formation in the rock-paper-scissors game 
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We investigate the impact of cyclic competition on pattern formation in the rock-paper-scissors game. By 
separately considering random and prepared initial conditions, we observe a critical influence of the competition 
rate p on the stability of spiral waves and on the emergence of biodiversity. In particular, while increasing 
values of p promote biodiversity, they may act detrimental on spatial pattern formation. For random initial 
conditions, we observe a phase transition from biodiversity to an absorbing phase, whereby the critical value 
of mobility grows linearly with increasing values of p on a log-log scale, but then saturates as p becomes 
large. For prepared initial conditions, we observe the formation of single-armed spirals, but only for values 
of p that are below a critical value. Once above, the spirals break up and form disordered spatial structures, 
mainly because of the percolation of vacant sites. Thus, there exists a critical value of the competition rate p c 
for stable single-armed spirals in finite populations. Importantly though, p c increases with increasing system 
size, because noise reinforces the disintegration of ordered patterns. In addition, we also find that p c increases 
with the mobility. These phenomena are reproduced by a deterministic model that is based on nonlinear partial 
differential equations. Our findings indicate that competition is vital for the sustenance of biodiversity and 
emergence of pattern formation in ecosystems governed by cyclical interactions. 

PACS numbers: 87.23.Cc, 87.23.Ge, 89.75.Fb 



I. INTRODUCTION 



Biodiversity and spatiotemporal dynamics of interacting in- 
dividuals are important for the characterization of ecologi- 
cal systems. Recently, spatial heterogeneity of species has 
attracted much attention because it is closely related to the 
stability and coexistence in ecological and evolutionary sys- 
tems |0]-|a]» whereas the competition between species was 
for example invoked as an evolutionary force in community 
food webs 10, 0]. Since each species has a different compe- 
tition rate with which it is able to invade others, Huntley and 
Kowaleski studied Phanerozoic fossil records and found that 
the competition rate is deeply interrelated with biodiversity 
JH]. Theoretically, several aspects of multi-species cyclical 
dominance have already been studied in detail. For example, 
it has been established that three species in a cyclic dominance 
exhibit self-organizing behavior on the spatial grid (HI fTTIl . 
whereby similar observations can be made also for system 
that incorporate more than three species, provided their to- 
tal number does not exceed fourteen [12]. Phase transitions 
and selection have also been studied in the predator-prey mod- 
els allowing some motion throughout inhabitable vacant sites 
El El. Recent theoretical studies focused on biodiversity 
and the emergence of self-organizing patterns by competing 
associations with homogeneous and heterogeneous invasion 
rates 1115143411 , elevating models of cyclical interactions to a 
very vibrant and enthusiastically approached field of research. 

Despite being fascinatingly simple and seemingly trivial, 
cyclical interactions among competing species emerge rather 
frequently in nature, and are indeed becoming more and 
more established as excellent models for the study of bio- 
diversity, the formation of defensive alliances [ 35] and Dar- 
winian selection 113611 . as well as structural complexity 113 711 



and prebiotic evolution in general [ 38] . Examples of re- 
ported real-life occurrences of cyclical interactions include 
rodents in the hight- Arctic tundra in Greenland IK39I1 . lizards 
in the inner Coast Range of California (400, overgrowths by 
marine sessile organisms lEHl . and microbial populations of 
Escherichia coll H42L l43ll . Recent experiments revealed that 
cyclical interactions can promote biodiversity within three 
strains of Escherichia coli JH). However, while effects of 
mobility and noise on cooperation, biodiversity and spatial 
pattern formation have been investigated intensely and with 
great success the effects of different rates of 

cyclic competition have not yet been thoroughly investigated. 
It is thus of interest to systemically investigate the impact of 
cyclic competition rate on both, biodiversity and pattern for- 
mation in models of cyclically competing species. 

Accordingly, we investigate effects of different competition 
rates on phase transitions and pattern formation by employing 
both random and prepared initial conditions. With increasing 
of mobility, the phase transition from biodiversity to an ab- 
sorbing phase, where two species go extinct, emerges, and the 
critical value of mobility is critically affected by the compe- 
tition rate. Large competition rates lead to vacant sites sur- 
rounding small patches containing all three species, which re- 
sults in the fragmentation of macroscopical spirals. This effect 
prevents the pattern outgrowing the system size, and in turn 
inducing promotion of biodiversity. For random initial con- 
ditions our results are in agrement with the seminal results of 
Reichenbach et al. DiSl, |T9|] . For prepared initial conditions, 
however, we find that with increasing of the competition ef- 
fects, a sharp transition occurs from a single spiral to multiple 
spirals coexisting. The disintegration is due to the percolation 
of vacant sites as well as noise. Since the level of the latter 
is directly related to the system size, the critical competition 
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rate p c increases with the increasing system size L x L, such 
that it is convenient to define the so-called competition effect 
as E = p/L, of which the critical value E c (at which the 
single-arm spiral disintegrate) is then independent of L. For a 
fixed system size, however, we also find that the value of the 
critical competition increases with increasing mobility. It is 
thus demonstrated that the rate of cyclic competition may have 
positive effects on the promotion of biodiversity, but negative 
effects on the sustenance of single-armed spiral waves. The 
latter phenomenon is investigated in detail also for different 
reproduction rates governing the rock-paper-scissors game. 

The remainder of this paper is organized as follows. In 
the next section we present the employed rock-paper-scissors 
game. Main results on effects of competition on biodiversity 
and pattern formation are presented in section 3. Section 4 
features results on the impact of different reproduction rates 
on the emergence and stability of spirals. The last section 
summarizes the findings and discusses their potential implica- 
tions. 



II. THE ROCK-PAPER-SCISSORS GAME 




FIG. 1 : (color online) (a) Extinction probability P ex t as a function 
of mobility M for different competition rates p. As M increases, 
there is a transition from stable coexistence (P ex t = 0) to extinction 
(Pext = 1). Panels (b) and (c) depict typical snapshots after a long 
relaxation period of the system for p = 1.0 and p = 10.0, respec- 
tively, at M = 1.0 x 10 -4 . (d) Phase diagram depicting the criti- 
cal mobility M c as a function of the competition rate p, separating 
the absorbing single- species phase (uniformity) and the biodiversity 
phase. Random initial conditions, N = 128 2 , and q = 1.0 were 
used for results in all panels. 



We employ the biological rock-paper-scissors game with 
the following specifications. Nodes of a L x L square lat- 
tice present mobile individuals belonging to one of the three 
species, which we denote by A, B and C. Each node can ei- 
ther host one individual of a given species or it can be vacant. 
Vacant sites, which we denote by <g>, are also the so-called 
resource sites. 
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where denotes any species or vacant sites. These reactions 
describe three processes, i.e. competition, reproduction and 
exchange, occurring only between neighboring nodes. In re- 
action (1), species A eliminates species B at a rate p, whereby 
the node previously hosting species B becomes vacant. In the 
same manner species B can eliminate species C, and species 
C can eliminate species A, thus forming a closed loop of dom- 
inance between them. Reaction (2) shows that individuals can 
place an offspring to a neighboring vacant node <g> at a rate q. 
Reaction (3) defines an exchange process, whereat an individ- 
ual may exchange its position with an individual belonging to 
any other species or an empty site at a rate 7. 

Here we use both the stochastic and the deterministic ap- 
proach to simulate the rock-paper-scissors game. For the 
stochastic approach, we make use of a stochastic simula- 
tion algorithm whereby the temporal evolution can be consid- 
ered as a random walk process. The most commonly applied 
stochastic simulation algorithm was developed by Gillespie 
lEH HHL where reactions occur in a random manner. In par- 
ticular, competition occurs with probability p/(p + q + 7), 
reproduction with probability l/(p + q + 7), and exchange 
(moving) with probability 7/(^ + ^ + 7). According to the 
random walk theory 115 3D , the mobility of individuals M can 



be defined as M = 7/27V, meaning it is proportional to the 
typical area explored by a mobile individual per unit time. 

As derived in the works of Reichenbach et al. fl8[[l9l l54ll. 
for the deterministic approach we use partial differential equa- 
tions (PDE) of the form: 

d t a(r, t) = DV 2 a(r, t) + qa(r, t)p - pc(r, t)a(r, t), 
dtb(r, t) = DV 2 b(r, t) + qb(r, t)p - pa(r, t)b(r, t), (4) 
dtc(r, t) = L>V 2 c(r, t) + qc(r, t)p - pb(r, *)c(r, t), 

where D denotes the diffusion rate and po the density of va- 
cant sites. 



III. EFFECTS OF COMPETITION ON BIODIVERSITY 
AND PATTERN FORMATION 

To examine effects of the competition rate p on the sus- 
tenance of biodiversity, we fix the reproduction rate q equal 
to 1.0 (see next section for the relaxation of this condition), 
and we perform extensive computer simulations of reactions 
(l)-(3) with no-flux boundary conditions, using first random 
initial configurations where each lattice site is occupied by an 
individual of species A, species B, species C, or is left empty 
with equal probability. Following Ref. Ill8n . we calculate the 
extinction probability P ext that two species have gone extinct 
when the system reaches the stationary state. Figure [TJa) fea- 
tures the results for p = 1.0, which are in agrement to those 
reported in Ref. Ill8ll . Namely, with increasing of M, a phase 
transition from biodiversity (P ex t = 0) to uniformity (P ex t = 
1) emerges at a critical value M c = (4.5 ± 0.5) x 10 -4 . From 
Fig -Ola), one can also observe that the critical value of M c de- 
pends onp; in particular, smaller p yield smaller values of M c . 
To study how the competition rate p affects the biodiversity, 
we present the typical snapshots for p = 1.0 and p = 10.0 
in Figs.QIb) and (c), respectively, at M = 1.0 x 10 -4 . We 
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FIG. 2: (color online) (a) Illustration of the realization of prepared 
initial conditions. Here R — 3.5, and every circle is occupied ex- 
clusively with a single species corresponding to the denoted color. 
All other sites are left vacant. Note that the distance between any 
two circles is the same. The radius R is fixed to 10.5 throughout 
this paper, (b) Phase transitions from single-arm spirals to coexisting 
multiple spirals, determined by means of P m , in dependence on E 
for M = 5.0 x 10 -5 , q = 1.0 and different system size N = L 2 . 
Here E describes the so-called competition effect, which we intro- 
duce as p/L. Prepared initial conditions, as depicted in panel (a), 
were used. 



find that while macroscopical spirals exist in the system at 
p = 1.0, at p = 10.0 small patches occupied by individuals 
of any of the three species are divided (or "disconnected") by 
vacant sites, and hence no macroscopical (large) spiral can be 
observed. According to Ref. [18], the loss of biodiversity re- 
sults from spirals outgrowing the system size when mobility 
M exceed M c . Therefore, for a given M, M = 1.0 x 10~ 4 
for example, the larger value of p promotes biodiversity bet- 
ter than smaller p. In Fig.QJd), we present the phase diagram 
with the biodiversity phase and the absorbing single-species 
(uniformity) phase delineated. There are in fact two differ- 
ent regimes inferable, which are due to small and large val- 
ues of p. For small values of p, i.e. when the reproduction 
rate is comparable to the competition rate, M c grows linearly 
with p on a log-log scale. However, for large values of p, i.e. 
when the competition rate is significantly larger than the re- 
production rate, the reproduction process limits the dynamics 
and therefore M c disobeys the linear relation and approaches 
a constant value for p > 10. 

Since pattern formation plays a central role in the suste- 
nance of biodiversity |[18l LL9D, we focus on the evolutionary 
process of pattern formation in dependence on the competi- 
tion rate p. Thus, we will consider p and the mobility M as 
the two crucial parameters effecting pattern formation in the 
examined model. As shown in Fig.[2ja), prepared initial con- 
dition that facilitate the emergence of spirals are employed 
for this purpose. In particular, three roundish areas with ra- 
dius of 10.5 are occupied individually by each of the three 
species, whereby the distance between each this region is the 
same. All other nodes are initially vacant, i.e. providing re- 
sources needed for reproduction. Notably, such a setup has 
been considered previously in ecological systems, for exam- 
ple those aiming for experimental bacteria growth 11551] . and is 
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FIG. 3: (color online) Panels in row (a) show single spirals at E — 
2.0 x 10 -3 for different system sizes. Panels in row (b) display 
spatial patterns at p = 2.5 for L = 200, L = 400, L = 800, and 
L = 1000. Accordingly, for L = 200 and L = 400, because the 
competition effect E is 1.25 x 10~ 2 and 6.25 x 10~ 3 , respectively 
(which exceeds E c ), multiple spirals are observed. However, for L = 
800 and L = 1000, since the competition effect E is 3.125 x 10" 3 
and 2.5 x 10 -3 , respectively (which is below E c ), single spirals are 
stable. Temporal evolution of typical spatial patterns, as emerging 
for p = 2.5 and N = 512 2 , is presented in row (c). At t = 100000, 
a single-arm spiral emerges when the system is started with prepared 
initial conditions. Subsequently, the large spiral starts to break up at 
t = 110900 and t = 111000 (indicated by the white arrow). When 
reaching the stationary state at t = 200000, the single-arm spiral 
observed at t = 100000 is fragmented. Here M = 5.0xl0" 5 and 
q — 1.0 in panels (a), (b) and (c). 



frequently referred to as growth initial conditions setup. 

We first investigate pattern formation by means of direct 
simulations of the game via the Monte Carlo method for dif- 
ferent values of p. An elementary Monte Carlo step consists 
of randomly choosing an individual who interacts with one of 
its four nearest neighbors, which is also selected randomly, 
and then executing the process as determined by the Gille- 
spie's algorithm |5j],[52|]. One full Mote Carlo step consists of 
N = L 2 elementary steps, during which, in accordance with 
the random sequential update, each player is selected once on 
average. Since noise increases with the decreasing of the sys- 
tem size, which on the other hand affects pattern formation, 
we define the so-called competition effect as E = p/L. In 
Fig. Ob), we plot the probability P m of multiple spirals coex- 
isting as obtained in dependence on E for M = 5.0 x 10 -5 . 
More accurately, 1 — P m is defined as the probability that the 
spatial grid is occupied by a single one-armed spiral. A sharp 
phase transition at a critical value E c = (4.5 ± 0.5) x 10 -3 
is observed. Because the transition becomes sharper with in- 
creasing system size, the phase transition appears to be dis- 
continuous. It is worth emphasizing that noise enhances the 
disintegration of spatial patterns, and thus the critical value 
of the competition rate p c increases with increasing system 
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FIG. 4: (color online) Typical spatial patterns emerging for dif- 
ferent competition rates p at M = 5.0 xlO -5 (top panels) and 
D — 5.0 x 10 -5 (bottom panels). Prepared initial conditions, as 
depicted in Fig. [2a) were used. From left to right, the values of p are 
0.1, 1.0 and 10.0, respectively. If p < p c = (2.3 ± 0.3) single-arm 
spirals emerge, which are characterized by decreasing spatial wave- 
lengths A = X/L as p increases. Upon surpassing the critical value, 
the large spirals disintegrate into several fragmented spirals, forming 
essentially disordered spatial portraits. The upper three panels are 
obtained from Monte Carlo simulations, while the lower three panels 
depict results from the PDE model. Here N — 512 2 and q = 1.0. 



size; for example, p c = (2.3 ± 0.1) for L = 512 while 
p c = (3.1 ± 0.1) for L = 700. Figure Oa) shows single 
spirals for different system sizes at^ = 2.0xl0 -3 [which is 
below the critical E c value that can be inferred from Fig. Ob)], 
while Fig. Ob) shows pattern formation for different system 
sizes at p = 2.5. Because competition effects at L = 200 
and L = 400 in Fig. Ob) yield E larger than E c , multiple 
spirals are observed, i.e., the large spiral disintegrates to sev- 
eral smaller and less ordered spirals. Competition effects at 
L = 800 and L = 1000 in Fig. Ob), on the other hand, 
are below E c , and therefore single-armed spirals are stable. 
Comparing the patterns at a given system size in Figs. Oa) 
and (b), one can also find that the wavelength of spirals de- 
creases with increasing of p. To get a better insight into the 
fragmentation of single-arm spirals, without loss of generality, 
we focus on the system size L = 512 where p c is (2.3 ± 0.3). 
In Fig.Oc), we first draw typical spatial patterns as obtained 
over time for p = 2.5 and M = 5.0 x 10 -5 . Starting from 
heterogeneous initial conditions, a single spiral emerges at 
t = 100000. However, one species in the arm of the single spi- 
ral breaks out at t = 110900 (marked by arrow), and shortly 
thereafter, two species break out at t = 111000. When the 
system reaches the stationary state at t = 200000, the single 
spiral is virtually completely fragmented and there is hardly 
any evidence left of its earlier existence. As shown in the top 
panels of Fig. ffl while globally ordered spiral waves emerge 
for p = 0.1 and 1.0, the latter disintegrate forp = 10.0. Fur- 
ther insights can be obtained by examining the spatial wave- 
length of spiral waves, defined as A = X/L, where X denotes 
the spatial distance between neighboring wave fronts with the 
same species. In particular, we find that the wavelength of 
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FIG. 5: (color online) Panel (a) shows the density of vacant sites po 
evolving with time for different competition rates p. Panel (b) shows 
the stationary values of po as function of the competition rate p, while 
the inset displays the data using log-log scale. Here M = 5.0 x 10 -5 , 
N 512 2 and q = 1.0. 



spiral waves decreases as p increases, and simultaneously the 
edges of the spirals separating different species become in- 
creasingly rough. We argue that these two facts eventually 
lead to the disintegration of globally ordered spiral waves for 
large enough values of p. Indeed, globally ordered spiral 
waves are no longer attainable for p > p c . This assertion 
is based on extensive Monte Carlo simulations revealing the 
presence of spatially periodic structures. This phenomenon is 
also reproduced by the PDE model [see Eq. (4)], as shown in 
the three bottom panels of Fig.[4j 

It is important to note that the competition rate p affects 
not only the density of the three species responsible for the 
emergence of spatial patterns, but also the density of vacant 
sites, as shown in Fig.Oa). With increasing values of p, the 
density of vacant sites increases as well, but the effect is of 
saturating nature, as can be inferred from Fig. Ob). For small 
values of p (at p = 10 -2 for example), the reproduction pro- 
cess (note that here q = 1.0) happens much faster than the 
competition process. There po grows linearly with p on a log- 
log scale (see inset of Fig. (b)). However, for larger p the 
reproduction process starts limiting the dynamics of po • Sim- 
ilarly as argued for results presented in Fig.QJd), po therefore 
starts disobeying the linear relation with p, and one can imag- 
ine that po will tend to 1.0 as p — » oo because the system will 
then stay in the prepared initial conditions. However, for very 
large values of p the density of vacant sites po rises slowly in- 
deed. In this case, the dynamics is almost entirely limited by 
reproduction, and patterns consisting of solitary wave fronts 
can be observed II 5 611. 

In Fig.Oa), we show further how the density of vacant sites 
varies with the critical competition rate p c as a function of mo- 
bility M. It can be observed that po remains small for all M 
at p = 0.01 because the reproduction process happens much 
faster than the competition process. However, po rises with 
increasing M for large p (at p = 10.0 for example), as there 
the reproduction process happens much slower than the com- 
petition process, and hence the large mobility M makes indi- 
viduals prey more effectively. In Fig. Ob) we finally plot the 
full p c — M phase diagram, where the biodiversity region is 
depicted white and the absorbing single-species phase (unifor- 
mity) is depicted shadowed. In addition, the line in the white 
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FIG. 7: (color online) Typical spatial patterns emerging for different 
reproduction rates g at M = 5. Ox 10 -5 . Prepared initial conditions, 
as depicted in Fig. 0a) were used. All patterns were observed after 
the system reached the stationary state. Here N — 512 2 andp = 1.0. 



FIG. 6: (color online) (a) The density of vacant sites po as func- 
tion of mobility M. (b) The critical competition rate p c as function 
of mobility M. The shadowed region denotes the emergence of an 
absorbing phase (uniformity), where two species go extinct. Here 
N = 512 2 and q = 1.0. 



region delineates single spirals (left) and multiple spirals co- 
existing (right). 



IV. EFFECTS OF REPRODUCTION RATE ON PATTERN 
FORMATION 

Based on thus far presented results, we can conclude that 
the cyclic competition has positive effects on biodiversity, 
similarly to effects of reproduction rate reported first by Re- 
ichenbach et ah 11811 . while it has negative effects on pattern 
formation. It is therefore of interest to test briefly whether the 
reproduction rate might also have negative effects on spatial 
pattern formation. To address this, we here fix the competi- 
tion rate to p — 1.0, and perform Monte Carlo simulations of 
reactions (l)-(3) for different values of the reproduction rate 
q. As above, we start from the prepared initial state depicted 
schematically in Fig. [2a), since it promotes pattern formation, 
in particular the emergence of spirals. 

Figure [7] presents typical spatial patterns emerging for dif- 
ferent reproduction rates q at M = 5.0 x 10 -5 . At q = 1.0, 
a globally stable spiral emerges, as we have already reported 
in the preceding section. With increasing values of of q , how- 
ever, the wavelength of single arm spirals first decreases (com- 
pare q = 1.0 and q = 10.0), and for higher q still (q=15.0 and 
20.0) the single arm spirals break up and become more and 
more fragmented. This is very much in agrement with what 
we have reported above for the impact of p by a fixed value 
of q, and it indeed confirms that increasing reproduction rates 
also negatively affect the emergence and stability of spatial 
patterns. 



V. DISCUSSION 

Based on a biological rock-paper-scissors game, effects of 
different rates of cyclic competition on spatial pattern forma- 
tion and biodiversity have been investigated. For low mobil- 
ity of individuals where the three species coexist, we have 



examined pattern formation affected by cyclic competition 
rates between species, discovering that self-organized struc- 
tures emerge if using prepared initial conditions. We have 
shown that globally ordered spiral waves, as were observed 
recently also in excitable systems |57|,[58|], can emerge from 
special heterogeneous initial conditions only if the competi- 
tion effect is smaller than a critical threshold E c . Since the 
introduced competition effect E = p/L depends on the sys- 
tem size, the critical competition rate p c is therefore itself 
system size dependent. This is because the system size de- 
termines the level of noise in the system, which in turn may 
facilitate the disintegration of spatial patterns. We found that 
when approaching p c at a given system size, the borders sepa- 
rating the species forming the spirals become more and more 
rough (non- smooth). This is because the density of vacant 
sites increases, but also due to noise that is inherent for finite 
system sizes. Accompanying this is also a decrease of spatial 
wavelength of spiral waves A, although the latter phenomenon 
bears no relevance for the impending disintegration of the spi- 
rals. At p = p c the disintegration of the globally ordered spi- 
rals occurs, resulting in predominantly disordered spatial por- 
traits consisting of small fragmented spiral-like patterns. For 
random initial condition, on the other hand, our results are in 
agrement with those reported in Ref. 1 18], where biodiversity 
affected by reproduction rate has been studied first. We have 
also investigated phase transitions D9I1 related to the extinc- 
tion process evoked by different competition rates, where we 
found that, similarly as by the reproduction rate, increasing 
cyclic competition rates might have a positive effect on biodi- 
versity due to the fragmentation of spirals, which prevents the 
patterns from outgrowing the system size. 

It is worth noting that spatiotemporal patterns have been in- 
vestigated extensively in the past in many different systems, 
ranging from chemical reactions on catalytic surfaces to prop- 
agating signals in aggregating microorganisms ll60ll . It has 
been shown that patterns in excitable systems emerge pri- 
marily due to the instabilities induced by the interplay be- 
tween the fast excitatory and slow recovery variables l6lll . 
This kind of mechanism explains well the spiral waves emerg- 
ing in the Belousov-Zhabotinsky (BZ ) reaction Ii62ll and ag- 
gregating amoeba D. discoideum 116311 . for example. Spirals 
in our systems, however, emerge because of the cyclic in- 
teraction between the three species, rather than differences 
in their dynamics. In addition, it is well known that the 
propagating signals and propagable interactions can lead to 
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complex spatiotemporal patterns in systems describing Ca 2+ 
signaling in thalamocortical neurons [64] or interactions be- 
tween predation and transport processes in a benthic nutrient- 
microorganism system [65]. Similar mechanisms can also re- 
sult in ordered spatiotemporal patterns occurring in the brain 
ll66l - [68ll and heart ll69ll . These mechanisms, however, are sig- 
nificantly different from the cyclic competition presented in 
this paper, and the observed patterns and waves are accord- 
ingly dissimilar too. For example, sequential waves were ob- 
served for brain Ii66l1 and heart 1 69] tissue, yet single-armed 
spirals as we report presently are rarely reported. Altogether, 
our findings thus indicate that the competition rate in models 
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